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ABSTRACT 

Observations of warped discs can give insight into the nature of angular momentum transport 
in accretion discs. Only a few objects are known to show strong periodicity on long timescales, 
but when such periodicity is present it is often attributed to precession of the accretion disc. 
The X-ray binary Hercules X-l/HZ Herculis (Her X-1) is one of the best examples of such 
periodicity and has been linked to disc precession since it was first observed. By using the 
current best-fitting models to Her X-1, which invoke precession driven by radiation warp¬ 
ing, I place a constraint on the effective viscosities that act in a warped disc. These effective 
viscosities almost certainly arise due to turbulence induced by the magneto-rotational instabil¬ 
ity. The constraints derived here are in agreement with analytical and numerical investigations 
into the nature of magneto-hydrodynamic disc turbulence, but at odds with some recent global 
simulations. 
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1 INTRODUCTION 

Accretion discs (Pringle & Rees 1972; Pringle 1981; Frank etal. 
2002) are found in many astrophysical phenomena, from star and 
planet formation to the gas accretion on to supermassive black 
holes which power active galactic nuclei. They also form in stel¬ 
lar mass binaries when, for example, a star fills its Roche lobe 
and loses mass to a neutron star or black hole companion. Her X-1 
is an X-ray binary in which a neutron star accretes from a donor 
star. This system displays some stochastic behaviour, but shows a 
strong 35-day periodicity (Tananbaum et al. 1972) linked to pre¬ 
cession of the accretion disc (Katz 1973; Roberts 1974; Petterson 
1975; Gerend & Boynton 1976). The precession is plausibly at¬ 
tributed to radiation warping of the disc (Pringle 1996; Pringle 
1997; Wijers & Pringle 1999; Ogilvie & Dubus 2001). 

Accretion discs evolve due to the action of a turbulent vis¬ 
cosity. In most cases this is driven by the magneto-rotational in¬ 
stability (MRI; Balbus & Hawley 1991). In some cases, e.g. mas¬ 
sive protostellar or AGN discs, the disc self-gravity can also in¬ 
duce turbulence alongside the MRI, but in X-ray binaries the disc 
mass is small enough that this does not contribute (e.g. eq. 5.51 of 
Frank et al. 2002). In the study of warped discs it is often assumed 
that the effective viscosity arising from the disc turbulence can be 
treated as isotropic. This means that azimuthal and vertical shear 
are damped at the same average rates. This effective viscosity is 
usually modelled as a Shakura & Sunyaev (1973) a viscosity. 

Warped discs are additionally complicated by the presence 
of a radial pressure gradient. This occurs because the midplanes 
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(and therefore regions of highest pressure) are misaligned if two 
neighbouring regions of the disc are misaligned. Therefore around 
each ring (It: in azimuth) there is an oscillating pressure force 
(Papaloizou & Pringle 1983). In the limit that a <K H/R this re¬ 
sults in a pressure warp wave being launched radially through 
the disc (Papaloizou & Lin 1995; Papaloizou & Terquem 1995; 
Lubow & Ogilvie 2000; Lubow et al. 2002). Ivanov & Papaloizou 
(2008) explore radiation warping dynamics for different disc mod¬ 
els, e.g. wavelike discs. If instead, a » H/R then the effect is 
strongly damped locally and the disc warp diffuses. X-ray binary 
discs are usually in the diffusive scenario as a ~ 0.1 (King et al. 
2007) and H/R ~ 0.01 (eq 5.50 Frank et al. 2002). 

In the diffusive warp propagation case, Papaloizou & Pringle 
(1983) provided the first self-consistent evolution equation for a 
warped disc, assuming an isotropic effective a. Their discovery of 
the oscillating radial pressure gradient showed that the effective vis¬ 
cosity governing the radial communication of the component of 
angular momentum parallel to the local orbital plane (v 2 = a 2 CsH) 
does not have the same magnitude as the effective viscosity govern¬ 
ing the radial communication of the component of angular momen¬ 
tum perpendicular to the plane of the disc (vi = a[CsH). In fact, 
while Vl oc Of, V 2 1 /a- 

The reason for this inverse dependence on the turbulent vis¬ 
cosity is the oscillating pressure gradient - the induced radial shear¬ 
ing motions are damped at a rate given by a and therefore commu¬ 
nication is restricted for larger values of a. A simple theory of this 
relation is given by Section 4.1 of Lodato & Pringle (2007), but 
the full derivation, in the linear and nonlinear cases respectively, is 
given by Papaloizou & Pringle (1983) and Ogilvie (1999). 

After the initial warped disc equations developed by 
Papaloizou & Pringle (1983), Pringle (1992) derived an evolution 
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equation for a warped disc based purely on conservation laws, 
which describes the evolution of a disc with an arbitrary shape 
once the effective viscosities v, and V 2 have been specified. Then, 
Ogilvie (1999, 2000) derived an evolution equation from first prin¬ 
ciples assuming little more than an isotropic effective a. Ogilvie’s 
equations are fully nonlinear and largely confirm the equations de¬ 
rived by Pringle (1992), but with two important differences: 1) an 
additional internal torque that makes tilted neighbouring rings pre- 
cess, 2 ) the effective viscosity coefficients ai and a 2 are functions 
of both a and the local disc warp amplitude. The work of Ogilvie 
(1999, 2000) gives a complete theory of the evolution of warped 
discs for an isotropic a in the diffusive case. For the wavelike case 
there is no full nonlinear theory (see e.g. Ogilvie 2006), but the 
linearised equations can be found in Papaloizou & Lin (1995) and 
Lubow & Ogilvie (2000). 

Almost all work on warped discs assumes that the effective 
viscosity arising from MHD turbulence is isotropic. Pringle (1992) 
noted that the viscosity may well be significantly anisotropic due 
to the different types of shear present in a warped disc. Azimuthal 
shear is secular where gas particles drift further apart, whereas ver¬ 
tical shear is oscillatory and thus should induce less dissipation. 
However, if the velocity spectrum of the turbulence is predomi¬ 
nantly on scales < H (e.g. Fig. 14 of Simon et al. 2012) then it is 
likely to act similarly in each direction. Theoretical investigations 
into the form of MHD disc turbulence are inconclusive, and so it is 
useful to examine what can be learnt from observations. King et al. 
(2013) argue that observations of warps can constrain the angular 
momentum transport in accretion discs, and in particular V 2 - This 
builds on King et al. (2007) who used observations of dwarf novae 
and soft X-ray transients to constrain the planar disc viscosity vj. 

In this work I explore the best-fitting models for the observed 
periodicity in Her X-1, which has a disc that is tilted and forced to 
process by the radiation warping instability, to see what constraints 
can be provided on the internal communication of angular momen¬ 
tum in a warped disc. 


2 WARPED DISC VISCOSITY 


The definitions and meanings attached to “viscosity” can be quite 
complex in the case of a warped disc, so I review the relevant 
equations here. The evolution equation for the angular momen¬ 
tum vector, L = 'LR^CU (where R is the radial coordinate, E 
is the disc surface density, O is the disc rotation law and / = 
(cos 7 sinjS, sin 7 sin/ 1 , cos/f) is the unit tilt vector for tilt and twist 
angles p and y respectively), in a warped disc (Pringle 1992; 
Ogilvie 1999) is 
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The first and second terms govern the diffusion of disc surface den¬ 
sity and tilt respectively and the third term is an advective term. 


The coefficients of the torques communicating angular momentum 
in the disc, vi and V 2 , are the effective viscosities. In reality these re¬ 
sult from MHD turbulence and hydrodynamics, but their effects are 
called viscosities, following a long-established usage in both the as- 
trophysical literature and the fluid dynamics literature. The second 
last term is the internal precession torque that is present in the anal¬ 
ysis of Papaloizou & Pringle (1983)', but first put into this form by 
Ogilvie (1999). In this term V 3 = a^CsH (where ~ 3/8; Ogilvie 
1999), but it is misleading to call this a ‘viscosity’ as it does not 
lead to diffusion, buf insfead causes a dispersive wavelike propaga¬ 
tion of the warp (Ogilvie 1999). As this term does not significantly 
affect the dynamics of diffusive warped discs, I do not discuss it 
further here. The final term corresponds to the external torque on 
the disc, in this case due to radiation. Pringle (1996, 1997) gives 
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where L* is the luminosity of the central source, c is the speed of 
light, the integral is taken only over illuminated parts of the disc, 
and A and are given by 


8p dy 

A = R — sm0 - R — cos 6 smyS 
8R ^ 8R ^ ^ 


and 


(3) 


= (cos cos 7 cos jS - sin sin 7 , (4) 

cos (j) sin 7 cos /I + sin cos 7 , 

- cos 0 sinjS). 

It can be seen from (1) that when there are no external torques 
on the disc (dG/dR = 0), the only steady disc shape is a flat disc. 
This occurs because the vectors /, dljdR and / x 81I8R form an 
orthogonal set. So the independent terms in (1) cannot balance each 
other, and so in the presence of dissipation both the second and 
fourth terms must decrease to zero, implying that \8l/8R\ = 0, i.e. a 
flat disc. 

Both effective viscosities, vi and y 2 , can be written in the 
Shakura-Sunyaev form 

vi = aiC^H; V2-a2CsH, (5) 

where ai and a 2 are the effective viscosity coefficients. It is often 
assumed that an “isotropic viscosity” is equivalent to taking these 
coefficients as equal, but this is not the case. Initial investigations 
into warped discs (e.g. Bardeen & Petterson 1975; Hatchett et al. 
1981) made this assumption, but it was shown to be inconsistent by 
Papaloizou & Pringle (1983). This is because of the internal oscil¬ 
lating pressure gradient described above. Instead, if we assume that 
a acts isotropically in the disc, damping any type of shear at the 
same rate, then Papaloizou & Pringle (1983) showed that 

1 

ai =a\ a2 = — . ( 6 ) 

2a 

This is what characterises an isotropic viscosity. Further, 
Papaloizou & Pringle (1983) give higher order terms in these coef¬ 
ficients and also Ogilvie (1999) extends this to the nonlinear warp 
regime. See Appendix A for a discussion. 

Now, if the disc viscosity is not isotropic then the notation 
becomes even more complex. We must now split (the previously 
assumed isotropic) a into two components tth and a^, which gov- 


* The imaginary part of the quantity A defined in their Section 4.2. 
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em the damping of azimuthaP and vertical shear respectively. For¬ 
mally tth corresponds to the horizontal-horizontal components of 
the rate-of-strain tensor, while o-y corresponds to the horizontal- 
vertical components of the rate-of-strain tensor (Torkelsson et al. 
2000). When cth = ffy = a (the isotropic case) we have effective 
viscosities given by (5) and (6), but when this is not the case we 
have 

O'! = Q\i\ 0^2 = ;;— • (7) 

ZQfv 

In principle these damping rates (orh and tty) can be calculated 
from MHD turbulence calculations and compared to determine if 
the effective a is isotropic or whether the azimuthal and vertical 
shear are not damped at the same rates. The initial investigation into 
this was the shearing box calculations of Torkelsson et al. (2000) 
that showed aj, » ay in good agreement with an isotropic a viscos¬ 
ity. This was further confirmed by the analytical treatment of MHD 
turbulence by Ogilvie (2003). 

In any case, it appears sensible to use observations to pro¬ 
vide any possible insight into the angular momentum transport in 
warped discs. In the next section I derive a constraint on the ratio 
of the vertical and azimuthal effective viscosities, by considering 
the best-fitting models to the tilt and precession in the X-ray binary 
system Her X-1. 


3 ANISOTROPY CONSTRAINT 

We consider the disc system in Hercules X-l/HZ Herculis (Her X- 
1). This is a semi-detached binary system with an F star of mass 
» 2Mq which transfers mass to a neutron star companion of mass 
~ IMg. The binary period is 1.7 days and the long period associated 
with precession is 35 days. Using appropriate parameters for the 
system. King et al. (2013) estimate the alignment timescale of the 
disc with the binary to be of order 10 days. Therefore, as the system 
has been observed to show this periodicity for > 40 years, it is clear 
that some physical process must be maintaining the tilt. 

Radiation warping (Petterson 1977b,c; Pringle 1996, 1997) 
can induce and maintain the tilt in such systems. This was first sug¬ 
gested for Her X-1 by Petterson (1977a) and further models have 
been developed by Wijers & Pringle (1999) and Ogilvie & Dubus 
(2001). Radiation warping is also successful in reproducing the 
long period properties of a number of other binary systems, as 
discussed in Wijers & Pringle (1999) and Ogilvie & Dubus (2001). 
The models have some unique and appealing properties consistent 
with a variety of observed effects, which I list here. The radiation 
warping mechanism is one of only a small number of physical ef¬ 
fects that can produce the disc tilt, and once the tilt is induced it also 
drives precession at a rate consistent with observations - in com¬ 
parison, precession forced by the companion’s gravity does not re¬ 
produce the observed precession rates (Wijers & Pringle 1999) and 
does not account for the disc tilt as the disc would align on a short 
timescale (King et al. 2013). The precession induced by radiation 
warping can be prograde or retrograde and therefore account for the 
unusual prograde precession observed in systems such as Cyg X- 
2 (e.g. Maloney & Begelman 1997; Brocksopp et al. 1999). If the 
central source luminosity is too high, the warped disc structure can 


^ The origin of this notation is from shearing box calculations which dis¬ 
cuss “horizontal” and vertical shear. In a full 3D disc, this corresponds to 
“azimuthal” and vertical shear. However, we are left with at, for azimuthal 
shear'. 


be non-stationary with oscillating tilt - this may be applicable to 
many of the long X-ray periods which are not stable in amplitude 
or period (Wijers & Pringle 1999). Another compelling feature of 
radiation warping is that the inner disc may tilt through more than 
90°, allowing it to accrete on to the central object in a counterro¬ 
tating manner - this has successfully explained the torque rever¬ 
sals seen in some neutron star X-ray binaries (van Kerkwijk et al. 
1998). These attributes of the radiation warping model form a com¬ 
pelling case for it being responsible for the long periods in a number 
of systems such as Her X-1. 

Using the models for Her X-1 provided by Wijers & Pringle 
(1999) and Ogilvie & Dubus (2001) we can constrain the effective 
viscosity (angular momentum communication) in the disc through 
the appropriate instability condition, as the disc must be unstable to 
radiation warping (but only just so) to undergo repeated precession 
(Ogilvie & Dubus 2001). Pringle (1996) gives the criterion for a 
disc to be unstable to radiation warping as 


R 


( 8 ) 


where R is the disc radius, R^ = GMtjd is the Schwarzchild 
radius of the accretor, zy = V 2 /V 1 is the ratio of vertical over 
azimuthal effective viscosities and e is the accretion efficiency. 
Ogilvie & Dubus (2001) write this criterion in terms of the binary 
orbit as 


Rb ^ 16 ;r-? 7 ^ Rb 

^ ~ 'R' 


For Her X-1, Rb/Rg = 3.1 x 10^, and the disc is expected to truncate 
at R » 0.3Rb for this near-equal mass binary with Mq/Mi = 1.56, 
where Mq is the companion mass (Ogilvie & Dubus 2001). Choos¬ 
ing R = 0.3Rb is a conservative limit as the instability must act at 
some radius in the disc R < 0.3Rb. For these numbers the criterion 
becomes 

-<75. (10) 

e 

Assuming a typical accretion efficiency e = 0.1 we now have a 
simple constraint on the ratio of effective viscosities. For radiation 
warping to drive the precession observed in Her X-1, we require 

Vb 

— <7.5. (11) 

Vl 

In the usual, isotropic, picture of accretion disc vis¬ 
cosity (Shakura & Sunyaev 1973; Papaloizou & Pringle 1983; 
Ogilvie 1999; Ogilvie 2000; Torkelsson et al. 2000; Ogilvie 2003; 
King et al. 2007; King et al. 2013) and in the limit of small tilt, the 
viscosity ratio can be written as zy k 1/2q'^ (Papaloizou & Pringle 
1983; see Ogilvie 1999 for the nonlinear corrections due to warp 
amplitude). In this case the limit is o- > 0.25, consistent with es¬ 
timates by Wijers & Pringle (1999), Ogilvie & Dubus (2001) and 
King et al. (2013). This limit is also consistent with the sugges¬ 
tion that Her X-1 is close to the stability limit (Ogilvie & Dubus 
2001), as ff w 0.25 (King et al. 2007)^. In other words, an isotropic 


^ Existing MHD simulations produce a values that are lower than inferred 
for Her X-1 and other similar systems. Shearing box MHD simulations find 
a ~ 0.01 (assuming no net vertical magnetic flux). For example Simon et al. 
(2012) find a few X 0.01. Global MHD simulations suggest numbers slightly 
higher, e.g. Fig. 1 of Parkin (2014) shows an n of ~ 0.05. However, there 
is currently no accepted reason for this discrepancy with observations (see 
e.g. King et al. 2007). The observations appear to be coirectly interpreted 
(e.g. Kotko & Lasota 2012) and the simulations appear to be converged (e.g. 


© 2014 RAS, MNRAS 000, 1-7 




4 Nixon 


(Navier-Stokes) viscosity is perfectly consistent with the idea that 
the precessing disc in Her X-1 arises from radiation warping (see 
also King et al. 2013). 

If we now instead assume that the effective viscosity is 
anisotropic, we have from (7) 

^ = ^ = ( 12 ) 

V] ffi 

(Torkelsson et al. 2000). Therefore the criterion for the radiation 
warping instability to act is 

:r^<7.5. (13) 

zorhQ'v 

If we again take ah(= «) = 0.25, then this becomes a 2 < 2, or 
equivalently tty > 0.25. Thus, for a significantly anisotropic vis¬ 
cosity, the action of radiation warping in Her X-1 places a limit on 
02 . 

There are higher-order corrections to this result, from both 
the limit of large a (Papaloizou & Pringle 1983; Kumar & Pringle 
1985) and the disc rotation law, which is modified by the binary 
companion. I discuss and derive the modification to this criterion 
due to both of these effects in Appendices A & B. I find that for the 
parameters in this study the corrections are relatively small, with 
less than a factor of two in the constraint. However, I note that for 
other systems these corrections can be large. 

Recently Sorathia et al. (2013) and Morales Teixeira et al. 
(2014) have argued that the effective viscosities in a warped disc 
arising from MHD turbulence are significantly anisotropic. For ex¬ 
ample, Sorathia et al. (2013) perform a single MHD simulation 
of a tilted accretion disc and report that tty =» 3 x 10“^ t ai,. 
Sorathia et al. (2013) do not report tth, but similar simulations yield 
values 0.01. Therefore the effective viscosity in Sorathia et al.’s 
simulation is significantly stronger than isotropy would predict. In 
their simulation this produces rapid disc alignment, on approxi¬ 
mately the (artificially shortened) precession time. Such disc be¬ 
haviour is not consistent with constraint (11) derived above from 
the precession observed in Her X-1. However, it is not clear that 
these simulations adequately resolve the MRI process (see e.g. 
Section 2.2 of Morales Teixeira et al. 2014). Also Sorathia et al. 
(2013) employ a precession timescale (almost as fast as dynam¬ 
ical) and disc parameters (H/R » a) that puts them into the 
highly nonlinear wavelike warps regime where there is no full 
nonlinear theory (even for an isotropic viscosity). Thus it is not 
clear how to interpret these results. Further, plots of components 
of the stress tensor as functions of position angles, as found in 
Sorathia et al. (2013) and Morales Teixeira et al. (2014) are some¬ 
what meaningless when comparing to an a theory, particularly 
whether the stress is positive or negative, as it is the average rate 
of damping that is given by a, and not the instantaneous local 
stress on scales <Si H x H. Thus it appears plausible that simply 
the wrong quantities are being measured from the simulations. On 
this point I finally note that Zhuravlev et al. (2014) report good 
agreement (in the wavelike warps case) between a semi-analytical 
approach based on an isotropic a and the MHD simulations of 
Morales Teixeira et al. (2014). They also report “minimal” differ¬ 
ences between an isotropic and anisotropic viscosity (see their Sec¬ 
tion 4.1). This presumably occurs because the discs are wavelike 


Simon et al. 2012; but see also Bodo et al. 2014). So the discrepancy seems 
real, and is therefore probably driven by either a missing piece of input 
physics or an environmental effect. 


and thus the short term dynamics are governed by pressure, rather 
than viscosity. 

Other investigations into the nature of MHD turbulence pro¬ 
duce results which are consistent with the criterion derived here. 
Torkelsson et al. (2000) measure tth and Oy from shearing box cal¬ 
culations and conclude they are both » 0.01, in agreement with the 
isotropic case. Ogilvie (2003) followed this up with an analytical 
investigation and found the same result. I note that the approach of 
Torkelsson et al. (2000) allowed a direct measure of the values of 
tth and Oy, thus offering a clean answer, without the need to disen¬ 
tangle the effects of global simulations with more complex physics. 


4 CONCLUSIONS 

I have derived a constraint on the anisotropy of the effective viscosi¬ 
ties in warped accretion discs. This constraint is found by consid¬ 
ering the best-fitting models for the periodicity observed in Her X- 
1. These models use the radiation warping mechanism of Pringle 
(1996). The constraint is satisfied by the widely used effective vis¬ 
cosities that result from a (near) isotropic a, but is in conflict with 
significantly anisotropic viscosities where the component commu¬ 
nicating misaligned angular momentum is significantly stronger 
than an isotropic model would imply. Recent MHD simulations 
have suggested the effective viscosity is significantly anisotropic in 
this manner, but this would imply that radiation warping in Her X-1 
is prohibited. Many other investigations into the nature of MHD tur¬ 
bulence have concluded that the viscosity is close to isotropic in the 
sense that azimuthal and vertical shear are damped at the same av¬ 
erage rates (Torkelsson et al. 2000; Ogilvie 2003; King et al. 2013). 
The constraint on the viscosity derived here (eqns 11 & 13) is con¬ 
sistent with this view. 
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APPENDIX A: HIGHER ORDER VISCOSITY 
COEEEICIENTS 


In this appendix I describe the relation between the coefficients a 2 
and Ui derived in Papaloizou & Pringle (1983), Kumar & Pringle 
(1985) and Ogilvie (1999). 

First, the (linear) derivation by Papaloizou & Pringle (1983) 
yields an equation for the disc tilt evolution (their 3.17) which is 
comparable to (1) above. Their equation contains the time deriva¬ 
tive and advection of the disc tilt with external forcing, equated 
with the diffusion terms (second and fourth terms of Eq. 1). They 
define the diffusive term through a complex quantity F given by 
(3.23) in Papaloizou & Pringle (1983). This quantity F with units 
of vSQ is defined by 


F = AR— (—) + fif—) , 
dR\RQl \Rnr 

A = vSO -I- (vIHji — iQ.)D , 

B = - id) + {v/Hl - iQ.)E, 

D = + inyiQ - v/H^) + 3ivniH^]IX, 


(Al) 

(A2) 

(A3) 

(A4) 


E = + vlH^)l{HlX) 


and 

X = lOv 


[hi"' H: 




(A5) 


(A 6 ) 


Here U^KRSi) is a dimensionless measure of the disc tilt and Hr 
and are defined by (3.21) and (3.22) of Papaloizou & Pringle 
(1983). As noted by Papaloizou & Pringle (1983) making the as¬ 
sumption that V is independent of the density p and that the 
gas is isothermal yields Hr = = H. With this and the 


Shakura & Sunyaev (1973) viscosity aSl = v/H-, we can simplify 
the above. In this case 

X = -t- 2ia ), (A7) 

which gives 
Ev 

E =- (A 8 ) 

a 

and therefore 


B = 0. 


(A9) 


Whereas 


D = Y.Q.H^ 


and so 


, —a^ -I- 3ia — 1 
Q;(f!! -I- 2i) 


A = EQ-H^fia ), 


(AlO) 


(All) 


where 


f(a) = 


6 ia + 2 + ija 
a + 2 i 


(A12) 


This is the form of the viscosity coefficient used by 
Kumar & Pringle (1985). The fact that B = 0 is a direct con¬ 
sequence of assuming isothermal gas. If the gas is not isothermal 
then B t 0 (see Kumar 1988). 

In the above equations a complex notation has been employed 
where the tilt is assumed to remain small (commensurate with the 
linearisation of the equations). In this case the unit tilt vector / = 
(tx, ly, h) is assumed to have L ~ 1 and then the equations can be 
solved for the complex tilt variable W = + ily It is therefore 

possible, with a suitable choice of coordinate system, to see that 
the real and imaginary parts of f{a) correspond to the coefficients 
of V 2 and V 3 , which act in the direction of dljdR and / x dljdR 
respectively. 

So if we simplify (A 12) we get real and imaginary parts given 
by 


Re[/(tt)] = 


1 4(1 -h 7 q'2) 
2 a 4 - 1 - 0 '^ 


and 


Im[/(Q')] = 


3(1 -2a^) 
4 -I- 


(A13) 

(A 14) 


Now If we note that it is the complex conjugate of A which en¬ 
ters the evolution equation (e.g. equation 2.1 of Kumar & Pringle 
1985) and that the version of the Pringle-Ogilvie equation (1) 
subsumes a factor of a half into V 3 , the relation between these 
coefficients and those derived by Ogilvie (1999) becomes clear. 
(A13) gives the relevant value of a 2 and -(A14)/2 gives the value 
of a^, as given by equation 137 of Ogilvie (1999) and the rela¬ 
tions Q '2 = 202 and a^ = Q^. Thus the a dependence derived 
by Papaloizou & Pringle (1983) and Ogilvie (1999) are identical. 
Ogilvie (1999) extends these coefficients to the nonlinear warp 
regime (see also Ogilvie 2000; Ogilvie & Latter 2013a,b), whereas 
Papaloizou & Pringle (1983) is valid for the linear regime. 

For small o- <K 1 we can see that these recover the classical 
values of a 2 and 03 relevant for the small warp and small a case. 
By ignoring terms of order a- we have 



and 



(A15) 


(Alb) 
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APPENDIX B: EFFECTS OF NON-KEPLERIAN 
POTENTIAL 


The disc in Her X-1 orbits in the potential of a binary system. 
Therefore the disc rotation, particularly in the outer regions is not 
strictly Keplerian. The warped disc evolution equation (1) accounts 
for an arbitrary rotation law Cl (Pringle 1992), and the radiation 
warping criterion is derived from this general form (Pringle 1996). 
The main effect of non-Keplerian rotation is to alter the values of 
the viscosity coefficients in (1). The instability of a disc to radia¬ 
tion warping requires an initially planar disc (with a small pertur¬ 
bation) to be linearly unstable (see e.g. Section 3 of Pringle 1996) 
and therefore it is appropriate to explore the linear warp viscosity 
coefficients. 

The leading order corrections to the coefficients, in the limit of 
small amplitude warps, is given by Equation A39 in Ogilvie (1999). 
This is 



i-2a + (1 - k^)ia^ 
(1 - k^) + 2ia - 


(Bl) 


where k = KjCl is the dimensionless epicyclic frequency and 
•240 = Qio + iQio is the viscosity coefficients evaluated at linear 
warp amplitudes (e.g. Q, = 2,o + |iAI^ Qa + f^di^l**)). By suitable 
rearrangement (set k = multiply by -1/ - 1, times by 2 and then 
take the complex conjugate) this recovers (A12). 

Now including these extra non-Keplerian terms, as above we 
have a 2 = 2Re[(24o]- So 


0-2 = Re 


i-2a + (1 - k^)ia^ 
(1 —k^) + 2ia - 


(B2) 


Simplifying and taking the real part gives 


2ak^ + 2(8 - k^)a^ 

(I -k^ - a~)^ + 


(B3) 


With ^ = 1 this recovers (A13). 

To evaluate this expression, for comparison with (A13) we 
must evaluate k^. The companion star induces perturbations to the 
disc orbits which deviate from Keplerian rotation. This effect is 
strongest at large disc radii closest to the companion. So for a con¬ 
servative estimate of the non-Keplerian effects we can focus on the 
disc outer edge at R = 0.3Rb- The equations for and 12^ in this 
potential are given by (e.g. Lubow & Ogilvie 2000) 


2 GM\ GM 2 
Cl - —^ -I- 


R3 2R^R 

and 

2 _ GM, GM 2 


Ai® 121 




1^1- 2b\‘’ I ^ 


( 1 ) 




where 


fo^”)(x)=- r 
^ ^ Jo 


cos(m0)(l + - 2xcos0) 


(B4) 


(B5) 


(B6) 


is the Laplace coefficient. 

Numerical evaluation at R = 0.3Rb, for an equal mass binary, 
yields k^ = 0.95. Evaluating (B3) with this value gives 0-2 = 2.78. 
This is close to the Keplerian value of 2.8 given by (A13). The dif¬ 
ferent calculations of q ;2 are shown in Fig. Bl This shows that the 
higher order a terms are more significant than the non-Keplerian 
terms for this system in determining a 2 . With both being close 
enough to the standard value that its use in Section 3 above is jus¬ 
tified. However, these additional calculations should be taken into 
account in more extreme systems, such as the inner parts of discs 
around Kerr black holes. 


3 


2.5 

s 

2 


1.5 


Figure Bl. Comparison of the three different calculations of the viscosity 
coefficient Q '2 for a = 0.25 and an equal-mass binary. The solid line is the 
standard (small a and Keplerian) value given by (6), the short-dashed line 
includes the higher order correction in a given by (A13), and the long- 
dashed line includes the non-Keplerian terms from the binary potential 
given by (B3). The assumed outer radius of the disc is shown by the ar¬ 
row on the x-axis at R = 0.3Rb. 
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APPENDIX C: HIGHER ORDER VISCOSITY 
COEFFICIENTS IN THE ANISOTROPIC CASE 


In Appendix AI discuss the higher order corrections to the viscos¬ 
ity coefficients assuming an isotropic viscosity. In this appendix I 
detail the same derivation, but retaining the anisotropic generality. 

Following Papaloizou & Pringle (1983) and noting that the R- 
0 stress coefficient is o-h, the R-z stress coefficient is ay and the (jt-z 
stress coefficient is (previously unintroduced) at (corresponding to 
Q '3 above), it can be seen that their (3.21) corresponds to ay and 
(3.22) corresponds to at. From this, and inspection of their (3.3) 
and (3.19), then (cf. Appendix A) 

A = VyZCl + (Vy I Hi - iCl)D (C1) 

where I have adopted the notation v, = ajH^Cl. Also, 

D = l.ClH\(vtlHl + iCl){iCl - Vb/H^) -b 3iVhn///^]/X (C2) 

and 

X = + ;^| + y.yJiHlHj). (C3) 

Assuming the same simplifications as before, this simplifies to 
X = UyatCl^ + iCl^{ay + art), (C4) 


D = 


[(o-t + i){i - oTh) + 3/o-h] 

OTvO't -H /(oTy -t- n't) 


(C5) 


and 


A = ILCl^H^ 


CYy, + ' 


(ay - ;■)([«, -I- /][; - orb] + 3iarh) 


O'vO't -I- i(ay + n't) 

This gives 

(oTv - -I- /][i - OTb] -b 3/Q'h) 


f(a) = ay + 


(C6) 


(C7) 


orvOr, -b i(ay -b ort) 

and therefore the real part gives 

Q'2 = (alal - atta^al + al + 3ayal + QrbQ,^ -b 2alatt (C8) 

-b 2ajat -b SQhQvQt + Qy + + a^ + af + 2ayQrt) . 
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A limit on disc viscosity anisotropy 


This can be understood by taking various limits. By setting 
ffh = ttv = cTt = a it can be seen that this recovers (A 13). If we 
assume all of the coefficients are small (dropping terms of order 
a]) and Ui ~ a^, this becomes 



(C9) 


Similarly, assuming all o-,- ~ 1, then 02 ~ 3.2. 

Note that in this anisotropic case, 6^0. This means that 
the diffusion of tilt governed by the V 2 term in (1) becomes an 
advection-diffusion term (see Kumar 1988) and this must be in¬ 
cluded for a complete picture. 
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